Excitonic absorption in gate controlled graphene quantum dots 
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We present a theory of excitonic processes in gate controlled graphene quantum dots. The de- 
pendence of the energy gap on shape, size and edge for graphene quantum dots with up to a 
million atoms is predicted. Using a combination of tight-binding, Hartree-Fock and configuration 
interaction methods, we show that triangular graphene quantum dots with zigzag edges exhibit 
optical transitions simultaneously in the THz, visible and UV spectral ranges, determined by strong 
electron-electron and excitonic interactions. The relationship between optical properties and finite 
magnetic moment and charge density controlled by an external gate is predicted. 
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Two-dimensional graphene monolayer exhibits fasci- 
nating electronic llHSj and optical properties |7H15| due to 
the zero energy gap and relativistic-like nature of quasi- 
particle dispersion close to the Fermi level. With recent 
improvements in nanofabrication techniques [16| the zero 
energy gap of bulk graphene can be opened via engi- 
neering size, shape, character of the edge and carrier 
density, and this in turn offers possibilities to simultane- 
ously control electronic [1y- 25], magnetic [161 1231430] and 
optical [30- 32] properties of a single- material nanostruc- 
ture. 

In this paper, we present a theory and results of numer- 
ical calculations predicting the dependence of the energy 
gap on shape, size and edge for graphene quantum dots 
with up to a million atoms. We show that triangular 
graphene quantum dots with zigzag edges combine mag- 
netic and optical properties tunable with carrier density, 
with optical transitions simultaneously in the THz, vis- 
ible and UV spectral ranges. We describe one electron 
properties of graphene quantum dots with N atoms and 
A^e ^z electrons by a combination of tight-binding ap- 
proach with a self-consistent Hartree-Fock method (TB- 
HF) described in detail in the EPAPS document of our 
earlier work Ref.(26). Then, in order to take into account 
correlation and excitonic effects, we solve the many-body 
Hamiltonian given by 



^ =Y^ V^pv V^ + Y^ ^pK<tK 



(1) 



p' ,(T 



p,a 



+ ^ E {p'l'\Vee\r's')bl,^bl,^,br'a'bs^c 

p' q r' s' 

+ 9 YkV(l\yee\rs)h\,Jl\,^,hr>a'hs'a 



pqrs 



Y. {{rp'\Vee\s'q) - 5,a'{rp'\VeMs'))bl,,h\^,hr,,bs,c 



p qrs 



where b^^,^ and /ij^ are hole and electron creation opera- 
tors corresponding to TB-HF quasi-particles. Excitonic 



100 



10 



LU 



0.1 



0.01 



'«•«•«•*•••«•• 



••••••••••••••• 



••••••••••••••• 



triangular • • •• 
zigzag 



iiexagonai 
armchair 



••••••••••• 

••••••••••••• 

••••••••••••••• 

••••••••••••••• 

••••••••••••• 

••••••••••• 

•••♦••••• 

-♦— Iiexagonai 
zigzag 




10° 10' 10' 10' 10' 10' 10' 

number of atoms 



FIG. 1: (Color online) Tight-binding bandgap energy as a 
function of total number of atoms N for a triangular zigzag 
quantum dot (black squares), hexagonal armchair quantum 
dot (red circles), and hexagonal zigzag quantum dot (blue di- 
amonds) . The inset shows the tight-binding energy spectrum 
in the vicinity of the Fermi level for the hexagonal armchair 
dot. 



absorption spectrum between ground state \ug) and ex- 
cited states Ivf) can then be calculated using 

A{oj) = ^ liuGlPWfMoJ - [Ef - Eg])) (2) 

/ 

where P = ^pp, Scja'{p\^\p^)hpabp'a' is the polarization 
operator. 

The electronic properties of graphene quantum dots 
depend on the size, shape and the character of the edge. 



This is illustrated by comparing electronic properties of 
three graphene quantum dots including (i) hexagonal dot 
with armchair edges, (ii) hexagonal dot with zigzag edges, 
and (iii) triangular dot with zigzag edges (see the inset 
of Fig. 1). The electronic structures are computed us- 
ing tight-binding Hamiltonian only with nearest-neighbor 
hopping for different number of atoms N. An exam- 
ple of the energy levels for A^ = 98 hexagonal quantum 
dot with armchair edges is shown in the inset of Fig. 
1. The red arrow indicates the bandgap separating the 
occupied valence band states from the empty conduc- 
tion band states. The dependence of the gap on the 
number of atoms is plotted in Fig.l. For the hexago- 
nal dot (red circles), the gap decays as the inverse of 
the square root of number of atoms N, from hundred to 
million atom nanostructures. This is expected for con- 
fined Dirac fermions with photon-like linear energy dis- 
persion {Egav^ kmin_r^ 27T / Ax oc l/vTV), as pointed 
out in Refs.|23|, Isjj, l33|). However, the replacement of 
the edge from hexagonal to the zigzag has a significant 
effect on the energy gap. The energy gap of hexagonal 
structure with zigzag edges decreases rapidly as the num- 
ber of atoms increases. This is due to the zigzag edges 
leading to localized states at the edge of the quantum 
dot, similar to whispering gallery modes of photons lo- 
calized at the edge of photonic microdisk [34]. Fig. 1 
also shows the effect on the energy gap of deforming the 
hexagonal structure into a triangle while keeping zigzag 
edges. The energy gap of a triangular dot is extracted 
from the energy spectrum, an example of which is shown 
in Fig. 2(a) for a triangular quantum dot with N=96. 
In addition to valence and conduction bands, the spec- 
trum shows a shell of degenerate levels at the Fermi level 
23|, l25|430| . The energy gap shown in Fig. 1 corresponds 
to transitions from the topmost valence to the lowest con- 
duction band state (see Fig. 2(a). Despite the presence 
of the zero energy shell, the energy gap in the triangular 
zigzag structure follows the power law Egap oc a/TV- We 
note that the energy gap changes from ^2.5 eV (green 
light) for a quantum dot with ^100 atoms to ^ 30 meV 
(8 THz) for a quantum dot with a million atoms and 
a diameter of ~ 100 nm. The presence of a partially 
occupied band of degenerate states in the middle of a 
well defined energy gap offers unique opportunity to si- 
multaneously control magnetic and optical properties of 
triangular graphene nanostructures. Indeed, due to the 
presence of the zero-energy band in the middle of the 
energy gap, several different photon energies correspond- 
ing to transitions within the zero-energy band, into and 
out of the zero-energy states, and valence-to-conduction 
band states are possible. This offers interesting possibili- 
ties for opto-electronic, opto-spintronic and intermediate- 
band solar cell photo- voltaic applications |32|, |35| . 

The dependence of optical and magnetic properties on 
the filling of the shell of degenerate zero-energy states 
is investigated in detail for a small triangular dot with 
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FIG. 2: (Color online) (a) Single particle tight-binding en- 
ergies of states near Fermi level for a N=96-atoni triangu- 
lar zigzag quantum dot. The colored arrows represent opti- 
cal transitions from valence to conduction band (VC, green), 
valence to zero-energy band (VZ, blue), zero-energy to con- 
duction band (ZC, red) and zero-energy to zero-energy band 
(ZZ, black). Inset shows corresponding dipole moments for all 
transitions up to 7 eV, obtained using tight-binding orbitals. 
Figs, (b-d) shows the effect of electron-electron interactions 
on the VZ transitions within (c) Hartree-Fock approximation, 
and including (d) correlations and excitonic effect obtained 
from exact configuration interaction calculations. 



zigzag edges containing N=96 atoms. For this system 
size, exact many-body calculations can be carried out. 
The degenerate band in the energy spectrum, shown in 
Fig. 2(a), has Nz = 7 zero-energy states. Each state 
is singly occupied and all electrons have parallel spin 
[26.-.29.] . We can thus classify allowed optical transitions 
into four classes, as shown in Fig. 2(a): (i) from va- 
lence band to zero-energy degenerate band (VZ transi- 
tions, blue color); (ii) from zero-energy band to conduc- 
tion band (ZC transitions, red color); (iii) from valence 
band to conduction band (VC transitions, green color); 
and finally, (iv) within zero-energy states (ZZ transitions, 
black color). As a consequence, there are three differ- 
ent photon energy scales involved in the absorption spec- 
trum. The corresponding joint optical density of states, 



calculated using dipole moments |(i|r|/)p connecting ini- 
tial and final states with energies ^^s and ^/s, are shown 
in the inset of Fig. 2(a). VC transitions (green) occur 
above full bandgap (^2.8 eV), VZ (blue) and ZC (red) 
transitions occur starting at half bandgap (~ 1.4 eV), 
and ZZ (black) transitions occur at THz energies. The 
energies corresponding to ZZ transitions are controlled by 
the second nearest-neighbor tunneling matrix element ^2 
and by electron-electron interactions. 

Fig. 2(b-d) illustrates in detail the effect of electron- 
electron and final-state (excitonic) interactions on the ab- 
sorption spectra. Fig. 2(b) shows detailed VZ absorption 
spectrum for noninteracting electrons. This spectrum 
corresponds to transitions from the filled valence band 
to half filled shell of N^ = 7 zero-energy states. Half- 
filling implies that each state of the zero-energy band 
is optically allowed. Numerical and analytical calcula- 
tions show that among the Nz = 7 zero-energy states 
there are two bulk like states, which couple strongly 
to the valence band resulting in the main transition at 
E = 1.41 eV. The remaining peaks in the absorption 
spectrum reflect transitions from the discrete spectrum 
of the valence band. When the electron-electron interac- 
tions are turned on, the ground state becomes a Slater 
determinant of self-consistently calculated Hartree-Fock 
orbitals: Valence states are doubly occupied, and zero- 
energy states are singly occupied by an electron with the 
same spin, e.g., down. The photon energies correspond- 
ing to optical transitions uj = {Ef + Sj) — {Ei + S^) 
are renormalized by the difference in quasi-particle self- 
energies Sj — S^. The absorption spectrum, shown in 
Fig. 2(c), is renormalized with transition energies blue- 
shifted by 0.51 eV to £; = 1.92 eV. Finally, when final 
state interactions between all interacting quasi-electron 
and quasi- hole states are taken into account, the excitonic 
spectrum is again renormalized from the quasi-particle 
spectrum, with transitions red shifted from quasi-particle 
transitions at ^ = 1.92 eV, down to E = 1.66 eV. As we 
can see, electron-electron interactions play an important 
role in determining energies and form of the absorption 
spectrum, with net blue shift from the noninteracting 
spectrum by 0.25 eV. 

We now turn to the analysis of the effect of carrier den- 
sity on the optical properties of graphene quantum dots. 
The finite carrier density, controlled by either metal- 
lic gate or via doping (intercalation), has been shown 
to significantly modify optical properties of graphene 
[Tl, Isl, UjJ, ll2| . For a quantum dot, the metallic gate shown 
in Fig. 3(a), changes the number of electrons in the de- 
generate shell from Nz to Nz + A A^^ . This is illustrated 
in Fig. 3(b) where 4 electrons were removed and 3 elec- 
trons remain. These remaining electrons populate degen- 
erate shell and their properties are entirely controlled by 
their interaction. Alternatively, removal of electrons from 
charge neutral shell corresponds to addition of holes. As 
is clear from Fig. 3(b), such a removal of electrons allows 
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FIG. 3: (Color online) (a) Schematic representation of a 
TGQD with N — 97 carbon atoms with 4 electrons moved 
to the metallic gate, (b) Corresponding single particle tight- 
binding configuration near the Fermi level, (c) Excitonic ab- 
sorption spectrum as a function of energy and charging ANz . 
For convenience, transitions are artificially broadened by 0.02 
eV. Peaks below 0.6 eV are due to ZZ transitions, peaks above 
1.2 eV are due to VZ (mostly on the left), and ZC (mostly 
on the right) transitions. Charge neutral case corresponds to 
ANz = (filling factors Ue = i^h = 1/2). At AN^ = 7, the de- 
generate band of zero states is completely filled with electrons 

(l^e = 1). 



intra-shell transitions ZZ, enhances VZ transitions by in- 
creasing the number of allowed final states and weakens 
the ZC transitions by decreasing the number of occupied 
initial states. Fig. 3(c) illustrates the overall effects in 
the computed excitonic absorption spectra for VZ, ZC, 
and ZZ transitions as a function of the number of addi- 
tional electrons AA/'^, related to shell filling and/or gate 
voltage. At AA^;^ = —7 (hole filling factor u^ = 1), the 
shell is empty and VZ transitions describe an exciton 
built of a hole in the valence band and an electron in 
the degenerate shell. The absorption spectrum has been 
described in Fig. 2(b-d) and is composed of one main ex- 
citonic peak at 1.66 eV. There are no ZC transitions and 
no ZZ transitions in the THz range. When we populate 
the shell with electrons, the VZ excitonic transition turns 
into a band of red shifted transitions corresponding to an 
exciton interacting with additional electrons, in analogy 
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FIG. 4: (Color online) Excitonic absorption spectrum at 
ANz =-1, 0, and 1 (on the left). Corresponding ground state 
spins are S = 3 (fully polarized), S = 7/2 (fully polarized), 
and S = (completely depolarized), respectively. The physics 
involved in optical transitions is schematically represented on 
the right side. 



to optical processes in the fractional quantum Hall ef- 
fect and charged semiconductor quantum dots [365 . As 
the shell filling increases, the number of available states 
decreases and the VZ transitions are quenched while ZC 
and ZZ transitions are enhanced. These results show 
that the absorption spectrum can be tuned by shell fill- 
ing, which can be experimentally controlled by applying 
a gate voltage. This is particularly true for the ZZ transi- 
tions in the THz range, which can be turned off by either 
emptying/filling the shell ANz = i7 or at half-filling. 
At half-filling, electron exchange leads to spin polariza- 
tion, with each state of the shell filled by a spin polar- 
ized electron. Since photons do not flip electron spin, no 
intra-shell transitions are allowed and the magnetic mo- 
ment of graphene quantum dot is directly reflected in the 
ZZ absorption spectrum. In addition, it was shown that 
even though at charge neutrality {AN^ = 0) the system 
is fully spin polarized, addition (but not subtraction) of a 
single electron leads to full spin depolarization [26]. This 
correlation effect plays an important role in the optical 
transitions involving zero energy states, shown in Fig. 
3(c) for ANz = ±1. 

In Fig. 4, we study the transitions for ANz = 0, ±1 
in detail. Fig. 4(b) shows the lack of absorption for 
half-filled spin polarized shell. The right hand side il- 
lustrates the fact that photons pass through since they 
are not able to induce electronic transitions and be ab- 
sorbed. For ANz = — I, Fig. 4(c), one electron is re- 
moved creating a hole in the spin polarized shell. Thus, 



the absorption spectrum corresponds to transitions from 
ground state to optically allowed excited states of the 
hole. The absorption spectrum for an additional elec- 
tron, ANz = +I7 shown in Fig. 4(a), is dramatically 
different. The additional electron depolarizes the spins 
of all electrons present, with total spin of the ground 
state S = 0[26]. The strongly correlated ground state 
has many configurations, which effectively allow for many 
transitions of the spin up and spin down electrons, as 
seen in Fig. 4(a). This asymmetry in the THz absorp- 
tion spectra allows for the optical detection of charge of 
the quantum dot and correlated electron states in the 
degenerate electronic shell. 
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